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Abstract 

We propose a new method for Hybrid Monte Carlo (HMC) simulations with odd numbers of dynamical fermions on 
the lattice. It employs a different approach from polynomial or rational HMC. In this method, 75 hermiticity of the 
lattice Dirac operators is crucial and it can be applied to Wilson, domain-wall, and overlap fermions. We compare 
HMC simulations with two degenerate flavors and (1 + 1) degenerate flavors using optimal domain-wall fermions. The 
ratio of the efficiency, (number of accepted trajectories) /(simulation time), is about 3:2. The relation between pseud- 
ofermion action of chirally symmetric lattice fermions in four-dimensional(overlap) and five-dimensional(domain- 
wall) representation are also analyzed. 

Keywords: Hybrid Monte Carlo, Odd Flavor Simulation, 75 Hermiticity, Schur Decomposition, Wilson Fermions, 
Domain-wall Fermions, Overlap Fermions 



1. Introduction 

In hybrid Monte Carlo simulations |[ll], the positive-definiteness of the action is essential to consider it as the 
statistical weight. When a lattice Dirac operator D is given, a positive-definite action of two degenerate flavors is 
easily constructed by using the hermitian conjugate of D, namely £)' D. The major difference of two flavor simulations 
and (2 -H 1) flavor simulations is that one cannot easily write down the pseudofermion action for the one flavor sector 
Rational or polynomial HMC methods , which approximates the square root of D^D, are mostly used for odd 
flavor simulations. 

In this paper, we provide a pseudofermion action for the one fl avor sector of the lattice fermions with 75 hermiticity 
without invoking the square root approximation for D. The main idea is very simple. For any lattice Dirac operators 
D with 75 hermiticity, P+DP+ and P_(1/D)P_ are hermitian, and one can construct a one-flavor pseudofermion 
action using these. The resultant action has the same determinant as D without any approximations. The non-trivial 
parts are to first check the positive-definiteness of the pseudofermion action and the discussion of how to obtain the 
pseudofermion action when there are mass preconditioners like the one in Hasenbusch method. 

The construction of the paper is as follows. In Sec.|2]-|5] the apphcation to Wilson fermions, Wilson fermions 
with mass preconditioner, domain-wall type fermions and overlap fermions are demonstrated, respectively. In Sec.|6] 
the relation between the pseudofermion action in four and five dimensional representations are presented. Numerical 
results are given in Sec.|2l while a summary and conclusion are provided in Sec.|8] 

2. Wilson Fermions 

In this section, we derive a pseudofermion action which will yield the same determinant as the Wilson-Dirac 
operator, 

r, I \ nu^ M_^V, ((W + m)l2x2 Yj^tf^cr^ \ 

Dwim) = (W + ot)1 -H > f^7^ = t \. , (1) 

V V tf,o-f, (W + m)l2x2 



where 
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l2x2 is two-by-two unit matrix and cr, (/ - 1,2, 3) are Pauli matrices. Throughout this work, we use chiral basis 
for gamma matrices. To get positive-definite pseudofermion action, we use a determinant relation which is valid for 
general matrix, 

det Dwim) ■ det[p'_Dwimy^ P'^) = det[p'^Dw{m)P'^) (2) 
where P'^i_ are the projectors which reduce two chiral components into one chiral sector, P'^ - (1 Q),P'_ = (0 1). 

The operators (P'_Dwimy^ P'_l) and (P'^Dw(m)P'_l) are hermitian due 75 symmetry. When the inverse of {W + m) is 
well-defined, the determinant of ([TJ is written as 



delDwim) - de.t[p'^Dw{m)P'l^ ■ det 



1 ^ 



P'_Dw{mT^P'l 



= det(W + mf det W//(m), (3) 



such that VK//(m) is the Schur complement of Dw{ni), i.e.. 



Wnim) = = (W -H m)l2x2 - V I^ttt^ — tyO-lcTy. (4) 

p'_Dw-^p1 7^^ + '" 

Thus, the pseudofermion action for one-flavor Wilson fermions can be written as 

Spp = (t|(W + mr^d), + (i>l(WHim))-^<i>2 = + m)-^<by - ^lP'_Hw{m)-^ P'l ^2 (5) 

where Hwim) - ysDwi'n), Oi is a pseudofermion field without a Dirac index, and <l>2 is a pseudofermion field with 
two spinor components. 

The positive-definiteness of the operator Wh('w) are discussed as follows. We note that for any background gauge 
field, the eigenvalues of W and (cr • t) satisfy the inequalitie^ < A(W) < 8, and \A(cr ■ t)\ < 4. It then follows that 
Wnim) is positive-definite for m > 4. Now, consider decreasing m from 4 to a smaller value. If the operator Wnim) is 
not positive for some m, e.g. m', then there must exist values of m at which det( W//(m)) is zero or singular in the region 
m' < m < 4. Among these values, if we denote the largest one as m„ , then for m > me,-, the positive-definiteness of 
WH(m) is assured. 

The value nicr corresponds to the opposite sign of the smallest eigenvalue of W. To see this, we use a relatioiJJ 

/l™n(W) < Re(A(Dwm) < A,^.AW). (6) 

From this relation, one can see that the value of m which makes det(W//(m)) singular (i.e. det(W -1- m) = 0) is 
larger than the value of m which satisfies det(Dw(m)) - det(WH(m)) = 0. There are no values of m which satisfy 
det (WHim)) = or det (WHim)) = +cx> above that. The condition for the positive-definiteness of VVh(w) can then be 
written as 

rricr — niQ s.t. det(Vt' + m*) — 0, mg > Wj > > ■ ■ ■ , 

0^ Wnim) (/) > for any (p and for m > nicr- (7) 

The smallest mass one can use in the one flavor method here is restricted by this bound. But it can be relaxed by using 
Hasenbusch's preconditioner as discussed in the next section. 

Generating Pseudo Fermion Fields 

To generate the pseudofermion field (^2 from a Gaussian random noise field H2, we need to take the square root of 
WH{ni), i.e., ^2 - V^h('w)S2 ■ This can be approximated by using the rational function of Wniin), 



O2 = /app (.Wnim)) 



Pi 



jrf ((qi + W + m)l2x2 - Z/i,v t^,^,tyO-^o-y ^ 



H2 - y/WH(m)E2, (8) 



Throughout this work, A(X) means any one of the eigenvalues of an operator X. 
^The proof of this relation is given in the appendix. 
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where po, pi and qi are expressed in terms of Jacobian elliptic functions. At first glance, the operations in ([8]l look 
formidable. However, since Wnim) is the Schur complement of Dwim), the inversion of the operator within the 
summation, in Eq. ^ can be obtained by the inversions of the operator (Dw(m) + qiP-), i.e.. 



+ ^ pi (Dw{m) + qiP^y 



1=1 



(9) 



Note that one cannot apply the multi-shift solver in Eq. (O, since P does not commute with Dw(m). However, for a 
given A^app number of inversions, the total number of iterations in the solver can be reduced by using the same idea as 
the chronological inversion method |4]. When one solves a set of linear equations {Dw{m) + qiP+)rii = H2 (1 < / < 
TVapp) with a given S2 serially from smaller I by using the iterative method, one can set a better initial guess rjf^ for the 
iterative method to solve {Dw{m) + P+)rii = H2 (2 < I) by using a linear combination of the solutions rjj (j < I) which 
have already been calculated i.e. rjf^ = Yjj<i <^iiVj- The coefficients cij are determined according to the prescription 
written in Ref. [(4:]. 

In this one flavor method, generating pseudofermions using the approximation given by Eq. (O makes the sim- 
ulation not exact. However, without using higher degrees of the approximation, one can make the algorithm exact 
by adding an accept/reject step after generating the pseudofemiion field. In the exact algorithm, the pseudofermion 
field should be produced according to a probability distribution proportional to e '^ wp^o"** xhis is obtained by mul- 
tiplying the operator a/WhC'w) to a Gaussian noise field H. In practice, however, the operator one uses in the sim- 
ulation is approximated by /appCWnC^M)) rather than VWh('w). This leads to a probability distribution that is pro- 

portional to e App(»'//i"'i)- . To adjust the difference, one can add an accept/reject step for with the probability 



/appiwH(m})2 / T\\K factor should be smaller than 1 and can be enforced by choosing /app(x) > ^]x for the 
whole eigenvalue region of Wnim). The discussion using eigenvectors is given in the appendix. 



3. Wilson fermions with the Hasenbusch method 

The idea of the Hasenbusch method [5} is to factorize the determinant det (Dw(nii)) into a product of determinants, 

detDwinii) - det{Dw(mi)/Dw(m2)) ■ det(Dw(m2)) (10) 

and the pseudofermion force coming from det (Dwimi) / Dw{m2)) is updated less frequently than the one coming from 
det(Dwim2)) in the molecular dynamics steps. The parameter 1112 is chosen such that the simulation cost is reduced. 
Furthermore, as already mentioned, in the one flavor method for Wilson fermions presented here, the factorization of 
the determinant allows us to use a smaller fermion mass in the simulation. For the second factor det{Dw(m2)), the 
same argument as the previous section Sec. ([2} can be applied. Here, we assume that mi < 1112 and that m2 is large 
enough such that the positive-definiteness of W//(m2) is always assured. 

Now, we consider how to treat the first factor det{Dw(mi)/ Dw(ni2)). Naively, one might consider the pseud- 
ofermion action like, 

-I- / W + OT2 \ , ,t (wH(m2) ^- \d>2 (11) 

^^\W + mij WH(mi)l 

but the second term is not hermitian in general, because W//(OTi)W//(m2) W//(m2)VV//(mi) for nii + m2- To remedy 
this situation, we rewrite the determinant d?X(D^(m\))j{D^ff(m2)') by using an operator which has different masses for 
the chkality plus and chirality minus sector, i.e.. 



detl^^ 

£>w('«2) 



detl (Dw'(O) H- P+mi -H P-m2) ^- | ■ dei\Dw{m2) ^- | - (12) 



Dwimi); \ DwiO) + P+rrii + P^m2 ) 



One may note that the first and the second factor on the right hand side are the same as the determinant of, 

l+(m2-m^)P'_-—^P'l \+{m2-m,)P' ^ -P'l (13) 

Dyf{m\) Dy^iS)) + miP+ + m2P- 
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respectively. Then, the pseudofermion action is written as 

SPF = 1>1<1'3 + <1>>4 + (m2 - mi)<t>lP' I P:'g>3 + (m2 - mi)<:>lP'^-—— ^- (14) 

^ '* J Dw{m\) DwiO) + m]P+ + m2P- 

where ^t, and 04 are pseudofermion fields with two spinor components. 

The condition for the positive-definiteness is given as follows. When nii - mi, the operators in Eg (fT3ll are 
the identity operators and their positive-definiteness is trivial. When one decreases m\ with OT2 fixed, the positive- 
definiteness will be lost only if either of the determinants in Eq. (fTZb becomes zero or singular. Thus, up to the largest 
m\ which satisfies det(Dvi/(0) + m\P+ + niiP-) - 0, the positive-definiteness of the operators in Eq (fTJt are assured. 
Note that this mass value is smaller than the limit given in Eq©. This can be seen by using the eigenvalue relation 
which is similar to Eq (|6]l, 

^nnniW + mj) < Re(i(Dw(0) + miP+ + rmP-)) < A^,AW + nii). (15) 

The largest m\ which yields det(Diy(0)H-OTiP+ +m2P-) = is larger than the largest m which yields det Dyfini) - 0. 
This can be understood from 

/mi + m2\ rriT — mi 
ysiDwiO) + miP+ + miP-) = Hw [ ^ J + ^ (16) 

and the properties of spectral flows of Hw{m) \\^. Moreover, for some gauge configurations, det{Dw(m)) may not 
be zero for any value of m. But there must exist mi which satisfies det{DwiO) + m\P+ + miP ) - for any gauge 
configurations. 



4. Domain-wall fermions 

The domain-wall type fermion operator |7] can be expressed as, 

£)dwf(w) = MDw{-ma){\ + cL(m)) + {I - L{m)) (17) 
with L(m) = P+L+im) + P^L^im) such that, 

, , , \<s<N,, l<s<N,, 

L+(m),,j' = <^ ^ , L^{m)s,s' = i ^ (18) 

\-mOs',N,7 l-OTd.5'_i, s-Ns, 

where m is the (bare) fermion mass, and mo G (0, 2) is a parameter called the "domain-wall height". The label s is the 
coordinate in the fifth dimension. Throughout this work, we assume that the number of sites in the fifth dimension Ns 
is even. The constant c and the diagonal matrix co - diag{wj) specify the type of domain-wall fermion. The operator 
Ddvitim) is the conventional domain-wall fermion for c = and cjs - 1 . It becomes an optimal domain-wall fermion 
when c = 1 and co/s are tuned such that maximal chiral symmetry is obtained. 

To obtain the pseudofermion action for one flavor, we modify the operator Ddwf(m) by multiplying 1/(1-1- cL) 
from the right. 

Hpf 

£)^„f(m) = £)dwf(m)/(l + cL) = coDw(-mo) + M(m, c) (19) 

with 

1 - L(m) 1 - L+(m) 1 - L_(m) 

M(m, c) = — - = i-— P+ + -— P_ = M+im, c)P+ + M^(m, c)P.. (20) 

1 + cL(m) 1 + cL+(m) 1 + cL4,m) 

In the following, we suppress the argument c of M^-im, c) for simplicity. Using the Schur decomposition, the determi- 
nant of the operator £)^^j (m) can be written as 

det£)'d„f(m) = det [oj{W - mo) + M^(m)f det^Wnim) = det [(x)(W - mo) + M^(m)f det nV//(m), (21) 



where 



(22) 



^VH(m) = R5 [co(W - mo) + M+(m)]l2x2 - t^- 



1 



(23) 



co(W - ma) + M-(m) 

Here is the reflection operator in the fifth dimension, (/Js) , ^' = ^s.a'j+.s'-i , which is introduced such that "Wnim) and 
'Wnim) are hermitian. For optimal domain-wall fermions, one still can choose a>/s which maintain maximal chiral 
symmetry and satisfies oj^^+i^s - <^s- After incorporating the contributions of the Pauli-Villars fields, the fermion 
determinant for domain-wall fermions becomes. 



det£)dwf('w) _ det [(d{W - mo) + M+(m)]^ de^Wnim) 
det£)dwf(l) det [a)(W - mo) + M+il)f det 

In principle, one can use a pseudofermion action like. 



(24) 



^(separate) ^ [^(^ _ + ^^(j)] 



1 



[co(W - mo) + M+(m)]' 



{(o{W - mo) + M+(l)] Oi-KDI 



1 



'Wuim) 



<1)2-h<D^'Wh(1)<1>3 
(25) 

But it is known that in the two-flavor simulation of domain-wall fermions, it is effective when a pseudofermion action 
uses a single set of pseudofermion field to estimate both the light fermion and Pauli-Villars terms [SJ. Then, we use 
the same we used in Hasenbusch method for Wilson fermions. 

Using the Schur decomposition of [£)dwf(l) - (A^+(l) - M+(m))P+'\, we obtain the relation 



det {u{W - mo) + M+{m)f ■ det [^h(ot) + A_(m)] = det [m(W - mo) + M41)f ■ det ['Wh(I) - A+(ot)] . 



where 

[A+(m)],,,, = [Rs (M+(l) - MAm))]s,.' , [A-(m)],,,, 
The properties of these matrices are given in the appendix. Using 

1 



det 



1 -I- A_(m)- 



•det 



1 -1- A+<m) 



[Rs (M_(l) - M_(m))],,,, . 
we can write the inverse of (l24t as 
1 



"WH{m)\ [ "WniD-A+im) 
and it can be used to construct the pseudofermion action. Using ( IB.6I) . we can simplify (l28l l to 

detyi ■ detS 

where. 



def 



1 +/(m,l,c)(v%). 



1 



"Wnim) 



(R5V)s' 



, def 



1 +g'im, l,c)(v^). 



1 



nV„(l)-A4m) 



Vs' 

s.s' J 



(26) 



(27) 



(28) 



(29) 



(30) 



The constant g' and vector v are given in Eq. (IB.6I 1 and (IB.4b . In the following, the arguments of g' are suppressed for 
simplicity (g' = g'(m, 1 , c)). Note that the five-dimensional matrix in Eq. ( l28T l is reduced to a four-dimensional matrix 
in this expression. Thus we can write the pseudofermion action for one-flavor domain-wall fermions as 



SpF = <D{^<Di H-<D^S02 

= <D|(Di -^'<l)[p:(y%) 

+ (Dj(D2+g'O^P;(v'), 



1 



ysRsD'd^fim) 
1 



r5-R5£''dwf(l)-A+(m)P+ 



(/?5v),'P:'Oi 

v.'P;'02 



(31) 



where and 02 are the pseudofermion fields (on the four-dimensional lattice) with two spinor components. 

Now we assert that the operators in (l29t are positive-definite for < m < 1 . At m = 1 , they are equal to the identity 
operator, and thus are positive-definite. As m is decreased, the operators in (|29) will cease to be positive -definite only 
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if either of the determinants in (|29) becomes zero or singular Using [W. ^M_(l)^=] - [W. -i=M+(m)-\=] = 0, and 
the fact that the eigenvalues of -^M_(l)-^ and -^M+(m)-^ have non-zero imaginary part^for < m < 1, we 
immediately see that (a>(W - nii)) + M_(l)) and (oj{W - mo) + M+{m)) cannot have a zero eigenvalue for < m < 1. 
Thus the operators in ( |29] l are well-defined for < ot < 1 . Furthermore, since (|28} is equal to the determinant of the 
four-dimensional Dirac operator with the approximation for the sign function of //kernel - J^Dwil -H (1 - c)Djv) the 
determinant cannot be zero and it follows that the operators in ( |29l ) are positive-definite for < m < 1 . 



Generating Pseudo Fermion Field 

We now discuss how to approximate the inverse square root of the operators ^ and S when one generates the 
pseudofermion fields Oi and <^2 from the Gaussian noise. Focusing on Eq. ( l35T l. we start from the inverse relation of 
five-dimensional operators, 



l+(J3-a)A^=- 



1 



Multiplying by 5 ^ from the left and S from the right, one obtains 



1 + (y6 - aVdiag(0, ■ ■ ■ ,0,1)5' 



1 +(a-;6)A+=- 



1 



"Wnim) + aA+ 



1 



1 + (a -/3)g'diag(0, ■ ■ ■ ,0, 1)5"' — 



1 



-s '^'m. 



Here, Aiss' = for s < s' . Then, the relation, Eq. (l34l i. also holds for the sub-block (s, s') - (Ns, Ns), 



1 



^H(m)+/3A^ 



1 



nV//(m) + aA^ 



One can obtain a similar equation for 'Wnim) and A_. 

For the square root of a general positive-definite operator A, the rational approximation can be used. 



(32) 
(33) 



(34) 



(35) 



In the case of y[, one has to calculate 



1 + g' v ^ R^ 



1 



/Wh(»j) 

Each term of the summation in Eq. ( |36] | is 

Pi 



1 



R^v\ H,. 



1 +qi 
Pi 



"WHim) + A_ 



1 + g' v^Rs 7 R5V 

l+qr ^VHim) + ^A^ 



l+qi (l+qd 
Pi 



Rsv 



1 +qi 



+ Pig' v'Rs 



1 



^H{m) + qiA 



-R-^v 



(36) 



(37) 

(38) 
(39) 
(40) 



^ Recall that M(m) for a conventional domain-wall fermion is a difference operator with anti-periodic boundary condition. 
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Then, to obtain (pi one needs to calculate 

" 1 

f^Hi ~^oHi + y p;^'v+7?5— — — -r-^5vEi (41) 
with, 

1 + ^; (1 + qiY 1 + ^; 

The same method can be used for the operator S. 

5. Overlap Fermions 

In this section, we construct the pseudofermion action for one flavor overlap fermions ||9|, 

Dov(m) = 1 + m + (1 - m)ys&ign{Hw(-mo)). (43) 
This operator satisfies Ginsparg-Wilson Relation(GWR) ifioll . 

Dov(m)^Dov(m) = (1 - m^) (Dov(O)' + 0^0)) + = 2(1 - m^) (f +Dov(0)P+ + P-Dov(O)P-) + 4toI (44) 

The overlap fermion also possesses 75 hermiticity, and so the application is similar to Wilson fermions. We begin by 
breaking this operator into its chiral components. 

i^ov(m) - yp, j^^^^^^pr, P'_DMP'l) ~ \D-.(m) D^.{m)) ^^^^ 
The determinant is written as, 

det(Dov(m)) = det(D++(m))det(D__(m) - D_+(m)— i— D+_(m) ) (46) 

By using GWR, one can show that the operators on the right hand side are positive-definite for < m < 1, 
provided that Dov{m) itself is well defined. The pseudofermion action is written as 

^ ',^+c(,lp'_—L-P'^4,2 (47) 



Practically, one has to use a reflection/refraction 111 111 or topology fixing term ill2ll to treat or avoid singularities in 
Dov('«) related to the topological change. 

The authors of UJJ proposed a simulation method for one-flavor with overlap fermions. We now highlight the 



differences between their work and our work. In Ref. Ill3ll . the authors use GWR, and factorize detZ)ov('w) £'ov('m) 
into two parts. 

detDov(m)"'Dov(m) = detP;Dov(m)"''Dov(m)f ■ detP'_DWn)^' DWn)P'J (48) 

The difference of the first factor and second factor of the right-hand side comes from the topological zero-mode of 
Dov(m). Then, the one-flavor determinant is written as 

detDov(m) = detP;Dov(m)"''Dov(m)f ■ m'^^^-^-^ = detP'_Doy(myDo^(m)P'J ■ m^^--^*K (49) 

Here, A^+/- are the numbers of topological zero-modes with a definite chirality. By using GWR, one obtains the 
relation between these operators and Schur complement. 



P'_DUmfDUm)P'J = Z)__(m) \D^Jm) - D^+(m)--^—D+Jm)] . 



(50) 



In other words, in Eq. (l4&t . the determinant is factorized as. 



det(Dov(m)) - det(D++(m)) ■ det 



'P'_DUmy'Do,{m)P'_ 



D..(m) 



- det(D++(m)) ■ det 



2(1 + m)D-4m) - 4m\ 
D__(m) j 



(51) 



Later, we show that the second factor corresponds to the term in Eq. ( ISTT i. Due to the cancellation between the 
numerator and the denominator, the force for the ^ term in HMC is smaller than that of the S term. 

When one performs a HMC simulation with a topological fixing term, it is apparent that using the pseudofermion 
action with the factorization in Eq. ( |49] | is more effective than the factorization in Eq. ( BTT i. For the five-dimensional 
representation Eq. (ISTT i. using only the S term and choosing the fermion mass parameter m' to satisfy (1 - m') - 
(1 - m^)/m, one can perform the HMC simulation respecting the factorization in Eq. ( |49] l. 



6. The relation of the pseudofermion action with four and five dimensional representations 

The ratio of determinants, Eq. (l24l . is equivalent to the determinant of the effective four-dimensional operator, 

D4d(m) = i(l +m)+ ^(1 - m)75/(Hkernel(-'«0)), (52) 

where the function f(x) is polar or a rational function which approximates the sign function. The form of the function 
/ is determined by cj. The operator i/kemei is defined as, 

H^emAi-mo) = ysT. — T, rrz—, r- (53) 

2 - (1 - c)Dw{-m(i) 

By using the Schur decomposition, the one flavor pseudofermion action is written as 

^\K—^P'l<^i+% ^ ^02. (54) 

Examining two cases, c - Q and c - 1, we show that the five-dimensional expression Eq. (ISTT l reproduces Eq. ( l54l l. 
This equivalence is not only of theoretical value, but aids in practical simulations. When one generates 4>\ and (pz, from 
Eq. (ISTT l. one has to know the eigenvalue spectrum of and S, a priori, in order to apply the rational approximation 
for the square root function. However, in the four-dimensional case, the spectrum is already known. 

case 7.- c = 

In this case, the propagator of the five-dimensional fermion at the boundary s - \ or Ns yields the propagator of 
the four-dimension fermion lll4ll . 

def (1 - ?M)£>4d 1 + r5/(^^kemel(-'M0)) 

Del, + m = — = — + m (55) 

^- = B—^ — -R5B\ (56) 

Dch + m Ddvfim) 

Here, B is defined as 

B, = P-6i,, + P+6n,,s (57) 

This relation holds only for c = 0. Note that this operator satisfies chiral symmetry Dd, 75 + 75 Dd, = in the limit 
Ns —> 00, in which f(x) becomes the sign function. From Eq. (ISST i and Eq. (|56] |. 



^ ■ " P^ + (l- m)P^ 



DAdim) 



1 



P . (58) 

s=\,s'=N, 



After incorporating the pseudofermion field, the left-hand side becomes the first term of Eq. (l54t . The right-hand side 
equals to the J[ term of Eq. OTT l by substituting g' - (\ - m) and v ' = (0, ■ ■ ■ , 0, 1). 
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m 


Ni 








Acceptance 


Acceptance/A^*i|,°''*'* 


0.019 


1 + 1 


75(1) 


260(1) 


345(1) 


0.88(3) 


2.6(l)x 10-" 


2 


0.6(1) 


239(2) 


240(2) 


0.90(3) 


3.8 (2) X 10-" 


0.038 


1 + 1 


54(1) 


125(1) 


179(1) 


0.90(3) 


5.0 (2) X 10-" 


2 


0.6(1) 


112(1) 


113(1) 


0.91(3) 


8.0 (3) X 10-" 



Table 1 : Comparison of HMC efficiency for the 2-flavor and (1 + l)-flavor QCD with optimal domain-wall quarks. The step size for the gauge field 
Aroaugc is 0.007(0.010) for m = 0.019(0.038). while the step size Arpp for (1 + l)-flavor pseudofermions is 0. 14(0.20) for m = 0.019(0.038), which 
is 4 times larger than that for the 2-flavor case. Here, A'j™', and ^tcr'"'' average CG iterations for one trajectory (for generating 

initial pseudofermion fields, molecular dynamics, and their sum respectively). 



Next, consider the S term. By using Eq. (l35T l. one obtains, 

1 



1 H- (1 - m)v' 



-1 



1 -(1 -m)p; 



1 



Drh + 1 



p': 



where we used the relation D^dim) - (Dch + m)/(Dch + 1). 



case 2: c - 1 



In this case, there are relations for and five-dimensional operators 111511 . 

[B£)dwf(m)-'£)dwf(l)B"'"] = £>4^(m)-', 
[B£)dwf(l)"' £)dwf(m)B^] = D4d(m). 



By using the relations , 



1 



1 - (1 - m)v' 



nvH(i) 



B(l + cLim))-' = — ^ (P+v"' + P-v'Rs), 



1 + c'^m 

Mm) (1 + cL(m))B'^ = g(m, 1, c) y[Tc{P+v + P-R^v), 
one can show that Eq. (|60] |. Eq. ( l6Tl l and projector reproduce the operators in Eq. (l3ll . 



(59) 



(60) 
(61) 

(62) 
(63) 



7. Numerical tests 

We compare the efficiency of the HMC simulation for two-flavor and (1 -i- l)-flavor QCD with domain-wall type 
fermion with c - \ and ojs - 1, on a 12-' x 24 x 16(A^s) lattice. For the gluon action, we use Iwasaki gauge action at 
P - 2.30. In the molecular dynamics, we use the Omelyan integrator |16], and the Sexton-Weingarten method lll7n . 
The pseudofermion action for the two-flavor simulation is the one with even-odd preconditioning which is described 
in Ref. 1 1 811 . The time step for the gauge field, (Arcauge), is the same for both two-flavor and (1 + l)-flavor cases, 
while the time step (Arpp) for the pseudofermion fields in the (1 + l)-flavor case is four times larger than that for the 
two-flavor case. The acceptance rate is roughly the same for both cases. We use conjugate gradient (CG) with mixed 
precision for the inversion of the quark matrix (with even-odd preconditioning). The length of each trajectory is set 
to two. After discarding 300 trajectories for thermalization, we accumulate 100 trajectories for the comparison of 
efficiency. Our results are given in Table{T] We see that the acceptance rate is almost the same for (1 + l)-flavor and 
two-flavor simulations. If the auto-correlation time is the same, then the efficiency of HMC can be estimated by the 
total acceptance divided by the CG iteration number, and the efficiency ratio for two-flavor and (1 -H l)-flavor is about 
3 : 2. 
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8. Concluding remarks 



In this work, we presented one-flavor method for HMC. The largest difference of the method presented here from 
RHMC and PHMC is that the pseudofermion action yields the one-flavor determinant without any approximations. 
For the overlap fermion, the difference from the method in Ref is that this can be used even if GWR is not exact. 
For the lattice fermions with 75 symmetry, it is always possible to construct real (hermitian) pseudofermion action, but 
one has to careful about the positive-definiteness, since it depends on the type of lattice fermions used. For chirally 
symmetric fermions like domain-wall/overlap fermions, positive-definiteness is assured in the entire mass parameter 
region which can be used in two-flavor simulations. On the other hand, for Wilson fermions, there exist a bound for the 
smallest mass where the positive-definiteness is assured. If the bound is larger than the mass value which is intended 
to use, one has to add Hasenbusch mass preconditioner then one can push the mass smaller Comparison between 
two-flavor and (I + l)-flavor using domain-wall type fermion HMC simulation shows that one can increase the step 
size of (1 -H l)-flavor simulation while keeping same acceptance ratio. The reason might be that for (1 + l)-simulation 
it is effectively the same as using four-dimensional operator and the force from the bulk mode is completely cancelled 
between the light fermion and Pauli- Villars field, while for the two-flavor pseudofermion action used in the comparison 
in this work, the cancellation was done only partly. The bottle neck of this method is generating pseudofermion field 
from Gaussian noise and it should be tuned and improved. Not only HMC itself, one can use these pseudofermion 
action for reweighting method, for example, to adjust the strange quark mass or to see the effects due to the difference 
of up down sea quark mass by using existing configurations. This approach for one-flavor simulation should be 
investigated further numerically and theoretically. 
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Appendix A. Tlie Eigenvalue Relation 

We prove a relation between "eigenvalues of the real part" and "the real part of the eigenvalues" for a general 
matrix M. 

i™„(Re[M]) < Re[A(M)] < A^,,(Re[M]), (A.l) 

def 

where A(M) can be any one of the eigenvalues of M, Re[M] is the real part of M, Re[M] = (M + M')/2, /linin(Re[M]) 
and /l„iin(Re[M]) are the smallest and the largest eigenvalue of Re[M], respectively. 

Proof 

The eigenvectors of Re[M] form an orthonormal basis, 

Re[M]cf>i ^ hi^i, (4>i,4>j)^6ij. (A.2) 

Here, ( , ) means the inner product. 

The eigenvectors i/r,- of M can be expressed in terms of the <pi's as, 

Mi//i^An//i, iJ^i^Yjc\Ui- (A.3) 

Here, we set (i/f,, i/^,) = 1 . This yields 2; |cJ'*P = 1 • The inner product of and Mi/f, gives the eigenvalue Ai, 

i, = (iA/,MiA,). (A.4) 

def 

Then, we divide M into the real part Re[M] and the imaginary part Im[M] = (M - M' )/2, and substitute these into 
the Eq. (lA4l l. 

Ai = ((Ai, Re[M] (A,) + Im[M] i/^d- (A.5) 
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After taking the real part of this equation, only the first term remains, 

Re(^,) = (iA/,Re[M] 1^,)- 
By substituting Eq. ( IA.3I ) into the right hand side of this equation, one obtains 

((A;, Re[M] (A,) = J] (cf )' c)'\a, Re[M] 0,) = ^ jcj'f /;,. 

This is equal to or larger than smallest eigenvalue of Re[M], 

2 |cff h, - ^„i„(Re[M]) = 2 |cff (/z, - ^™in(Re[M])) > 0. 

I I 

Thus, by using Eq. ( IA.6I ) to Eq. ( IA.8b . it is proven that Re[/l(M)] cannot be smaller than /tniin(Re[A/]), 

^™„(Re[M]) < Rs[A(M)]. 
The other inequality in Eq. (lA.lb . Re[/l(M)] < A^^^{Re[M]), can be proven in a similar way. 



(A.6) 



(A.7) 



(A.8) 



(A.9) 



Appendix B. Eigenvalues and Eigenvectors of Matrix \±{m) 

Here, we derive the eigenvectors and the eigenvalues of A+(m) defined in Eq. dZTl ). We work explicitly on A^^ = 4. 
After obtaining the answers for A^, - 4, to convert the answers for general Ns is straightforward. 
Using the matrix M+{m, c), i.e.. 



M+(m, c) = (1 - L+)(l + cL+y 



def 



the matrix A+(mi,m2) = Rs{M+{m2) - M_(mi)] is written as. 



A+(oti,OT2) = ^5 



1 + 



1 + c^m2 



1 c^m2 -c^m2 cm2 ^ 

—c 1 c^m2 -c^m2 
1 c^m2 
-c 1 ) 



—c 

^2 



1 + i 1 

1 + cL+ c ' 



i + i 



(B.l) 



1 + c'^niy 



1 




—c^nii 


cmi ' 


-c 


1 




— c^mi 




—c 


1 






c2 


-c 


1 J 



= g(mi,m2,c)Q(c) 
Here, g(mi,m2, c) and 2(c) are defined as. 



g{mi,m2,c) 



(1 + c)(ot2 - Wl) 
(1 +c4mi)(l +c4m2)' 





-c5 




















-c 


l-c3 


c2 


-c 


1. 



(B.2) 



The matrix A+(m) in EqdZTl) is given by A+{m) = A+(m, 1). 
Eigenvalues and eigenvectors of Q are. 



Qui^Q (i = 1,2, 3), e V = /If-v with = + c"* + c^ + 1, 









ro> 








c 




1 







Ml = 







c 




1 




.0. 




.0. 







1 



c 
-c 

I i; 



Define unitary matrix S as 



5 = Mj . 
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(B.3) 



(B.4) 



(B.5) 



Here, u[, u\ and u'^ are the vectors made by orthonormalization of u\, ui and M3. Then, A+(mi, mi) is written as, 

A+(mi,m2) ^ g(mi,m2,c) S dmg(0,0,0, Ac) = g'(mi,m2,c) S diag(0,0,0, 1) 5^ (B.6) 

with g'(mi,m2,c) — Acg{m\,m2,c). 

The similar relation for A_(oti, OT2) = /?5[M_(ot2) - M-{mi)\ is given by multiplying from left and right to this 
equation. 

Appendix C. Exactness - Point of View with Eigenvalues 

Consider the eigenvalues and eigenvectors of the operator Wh- 

WHi^i=hit^i (C.l) 
The Gaussian noise field H, and pseudo fermion field 0, are expressed as hnear combination of the eigenvectors (fr, . 

H = <^ = 2c,iA,- (C.2) 

The path integral of is written as the product of integral w.r.t the coefficient C/. 

When bi has a distribution e~''i''', one can make the distribution proportional to e"'^' by multiplying by y/hi, Ci = 
y/hibi. 

On the other hand, constructing c, = fs^ppih^bi, the distribution become e ' AppC-.*^ 
The difference can be adjusted by accept/reject step with the probability 



(C.4) 



To take this factor as a probabiUty, this must be less than one, i.e., fwp{hi) > y/hi. This condition is understood that 



g ' Apptfi/)^ ' has a slightly broader distribution than e ' V*; and the accept/reject step makes the distribution narrow 
by suppressing larger c, . 

Appendix D. Linear Algebra 

Here, we remind the readers some relations of Unear algebra. 
• Schur Decomposition 



A B\_l \ 0\lA \h A-^B\_h BD-^\lA- BD-^C 0\( 1 

C Z) ~ CA-i 10 D-CA-'B 1 ~ 1 D U)"'C 1 



(D.l) 



The Inversion of an Schur Complement 

The Determinant of an Operator with Projector P± 

det P'^MP'^ = detCP+M + P.) = detiMP+ + P.) (D.3) 

det(MP+ + P_) = det(M(P+ + M'^P^)) = det(M) det(P+ + M-^P.) = det(M) detiP'_M-^ P'^) (D.4) 
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